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In this paper, the localized surface modes in a defective multilayer structure has been investigated. 
It is shown that the defective multilayer structures can support two different kind of localized modes 
depending on the position and the thickness of the defect layer. One of these modes is localized at the 
interface between the multilayer structure and a homogeneous medium (the so-called surface mode) 
and the other one is localized at the defect layer (defect localized mode). We reveal that the presence 
of defect layer pushes the dispersion curve of surface modes to the lower or the upper edge of the 
photonic bandgap depending on the homogeneous medium is a left-handed or right-handed medium 
(e.g. vacuum), respectively. So, the existence region of the surface modes restricted. Moreover, the 
effect of defect on the energy flow velocity of the surface modes is discussed. 

PACS numbers: 42.70.Qs, 41.20.Jb 



I. INTRODUCTION 

The advent of photonic crystal (PC) materials [j], 
excited a lot of interest toward the existence of sur- 
face modes (SMs) at the interfaces of such materials 
[S0, IS Hi- These modes remain localized around the PC 
surface in a manner similar to surface plasmons on metal 
slabs Q. In real materials, PCs are always finite in size 
and, therefore, SMs can exist [1, H, H, Q. The existence 
of SMs can directly affect the performance and efficiency 
of PCs in applications. Thus, it is important to study 
the SMs in PCs. The first experimental observation of 
such SMs came by Robertson et al ||. In the experi- 
ment, the authors employed a standard attenuated total 
reflection (ATR) setup, widely used for surface-plasmon 
observations on metals [HI EH- The majority of the sub- 
sequent theoretical studies focused on the existence of 
such SMs in various PC structures. It was found that 
the frequency and in many cases even the mere existence 
of SMs are strongly influenced by the way the periodic 
PC is terminated 0, H|. Nevertheless, the initial acute 
interest for SMs had somewhat subsided until recently. 
The need to understand and engineer SMs came back to 
light, when these deemed to play a key role in newly dis- 
covered PC phenomena. In particular, it was found that 
SMs can strongly influence the subdiffraction focusing 
properties of PC based slab superlenses [H, EH EH EH- 
Moreover, coupling to such SMs in PC subwavelength- 
width waveguides leads to a highly directional exit beam 
[iH E3- Despite the intensive research on PC surface 
phenomena [H, 13, S @] , one aspect of the SM propaga- 
tion remains unexplored. When the periodicity of PC is 
broken by introducing a defect into a PC, a defect mode 
will appear inside the photonic bandgap (PBG) due to 
change of the interference behavior of light, whose prop- 
erties would be determined by the nature of the defect. 
The introduction of defect layers in one-dimensional PCs 
can create defect modes within the PBGs, just as defect 
layers in semiconductor superlattices may result in elec- 
tron defect states in the band gaps. A natural question 



arises: what is the effect of defect on the SMs? 

In this paper, we investigate the effect of defect on the 
SMs at the interface between a multilayer structure and 
a homogeneous medium (the interface of the structure). 
We show that two different kind of localized modes can 
be created in the defective multilayer structure depend- 
ing on the position and the thickness of defect layer. One 
of these modes is a SM and has a peak intensity at the 
interface of the structure. The other one is localized at 
the defect layer which we call it the defect localized mode 
(DLM). Furthermore, it is shown that the existence re- 
gions of the surface modes and defect localized modes 
depends on the position and thicknesses of the defect 
layer. As well, we show that the energy flow velocity of 
the surface modes has not been considerably affected due 
to the existence of the defect layer in the periodic multi- 
layer structure. However, the energy flow velocity of the 
defect localized modes is affected by the position and the 
thickness of the defect layer. In Sec. II, we introduce the 
model of the system under consideration. In Sec. Ill, the 
properties of localized modes are studied. Finally, Sec. 
IV concludes with brief comments. 



II. FORMALISM 

We wish to describe localized modes that form at 
the interface between a homogeneous medium of low re- 
fractive index, uq = ^/eo/^cb aud a defective multilayer 
structure with layers of refractive indices ri\ — ^Je\^i\ 
and ri2 = ^£2^2 and thicknesses d\ and da. Here we 
adopt the typical values used in reference fl8l |: e = — lj 
Mo = — 1) £\ = 4, /ii = 1, £2 = 2.25, /12 = 1, di = 1 cm 
and d2 = 1.65 cm. We suppose that the defect layer 
with thickness dd and refractive index n-y located after 
T], (i] — 0, 1, 2, ...) complete periods (d = d±+ cfe) of mul- 
tilayer structure, hence r\ describes the position of the 
defect layer in the structure. We choose a coordinate 
system in which the layers have normal vector along OZ 
and we consider the propagation of monochromatic TE- 
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polarized waves described by [19( 

E = E y (z)e y e i(kpx ~ UJt \ 



(1) 



with the electric field E in the y direction. Here, ui is the 
angular frequency, k — uj/c is the vacuum wavenumber, 
f3 = %■ and k x is the x-component of the wave-vector 
of modulus k(z) = kn(z). We look for stationary solu- 
tions propagating along the interface which satisfy the 
following scalar Hclmholtz-type equation 
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In the periodic structure, the waves are the Bloch modes 
E y (z) = ip(z)e lKbZ where K b is the Bloch wavenumber, 



and ip(z) is the Bloch function, which is periodic with the 
period of the photonic structure (see Ref. In the 

periodic structure the waves will be decaying provided Kb 
is complex; and this condition defines the spectral band 
gaps of an infinite multilayer structure. For the calcula- 
tion of the Bloch modes, we use the well-known transfer 
matrix method [20]. To find the localized modes, we take 
solutions of Eq. ^ in a homogeneous medium and the 
Bloch modes in the periodic structure and satisfy the 
conditions of continuity of the tangential components of 
the electric and magnetic fields the interface of the struc- 
ture [2l| . In this way, we can obtain the exact dispersion 
relation k = fc(/3) for TE-polarized localized modes by 
numerically solving the following dispersion condition: 
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where go — kyW 



and k\ — k^J n\ — (3 2 . Here, 



A, B are the elements of the transfer matrix of the 
multilayer structure, A is the eigenvalue of the trans- 
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III. NUMERICAL RESULTS AND 
DISCUSSIONS 



(SM) (see Fig. E](a)), while the mode 2 is localized at the 
defect layer (DLM) (see Fig. Efc)). 

Our investigations indicate that the existence region of 
the SMs and the DLMs strongly depends on the position 
of defect layer. To show this, we plotted the dispersion 
curves of the localized modes in the first spectral gap as 
a function of (3 for different defect position in Fig. [3j As 
one can see from Figs. E{a),(b), when the position of the 
defect layer approaches to the surface of the structure 
(small 77), the coupling between the incident wave and 



In the following, we want to discuss the effect of de- 
fect layer on the dispersion properties of the localized 
modes. In the used structure, we assumed that the ho- 
mogeneous medium is a left-handed metamaterial (LHM) 
with no = — 1. For comparison, we also investigated the 
dispersion of the corresponding localized modes in the 
structure, where the homogeneous medium is replaced by 
vacuum (no = 1). To do this, we plotted two sides of the 
dispersion condition (Eq. ([3])) vs. k in the first spectral 
bandgap in Fig. [T] In this figure the solid and dotted 
lines show the left-hand side of Eq. for the cases 
no = — 1 and no = 1 respectively, while the dashed lines 
represent the right-hand side of Eq. ([3]). Here we take the 
following values: dd = 0.8 di, r) = 3 and (3 = 1.14 cm -1 . 
Unlike the defectless structure we see that the dis- 
persion condition has two solutions (points 1 and 2 in 
Fig. [T|). So, there are two different modes for a given 
p. To describe these modes, we plotted the intensity dis- 
tribution of these modes as a function of z in Fig. [2] 
As one can see, both of these modes are localized modes. 
But, the mode 1 is localized the interface of the structure 
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FIG. 1: Graphical representation of the solutions of dispersion 
condition (Eq. (O). Here the solid and dotted lines show the 
left-hand side of Eq. ((3]) for the cases no = — 1 and no = 1 
respectively, while the dashed lines represent the right-hand 
side of Eq. ([3J. In our calculations we take the following 
values: dd = 0.8 di, f] = 3 and (3 = 1.14 cm -1 . 
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FIG. 2: Geometry and examples of the surface and defect 
modes: (a) SM, k = 0.821 cm -1 , (b) geometry and (c) DLM, 
k = 1.019 cm" 1 . These modes are corresponding to the points 
1 and 2 in Fig. [T] respectively 

. Here, (3 = 1.14 and no = — 1. 



the defect mode strongly take places and leads to the 
formation of the DLM. But by increasing n, the interac- 
tion between the incident wave and defect mode becomes 
weaker (see Figs. [3(c)), so that in the case of sufficiently 
big T) the DLM disappears in the benefit of the appear- 
ance of the SM. This case (77 corresponds to the 
defectless structure (see Figs. [3(d)). 

In Fig. [3J we also studied the dispersion property of 
two different homogeneous media (i.e. LHM and vac- 
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FIG. 3: The effect of defect position on the dispersion proper- 
ties of SMs (solid lines) and DLMs (dashed lines) in the first 
band gap. Here, the thin lines show the case for no = — 1 
and the thick lines show the case for no = 1. The shaded 
regions show propagation bands of the multilayer structure. 
Here, d d = 0.8 di. 
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FIG. 4: The calculated ATR spectrum of a defectless multi- 
layer structures with (a) (3 = 1.09, (b) f3 — 1.29 and a de- 
fective multilayer structures with (c) /3 = 1.09, (d) /3 = 1.29. 
Here, the thin lines show the case for no = — 1 and the thick 
lines show the case for no = 1. The other parameters are 
dd — 0.8 di, n = 3 and n = 3. 
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FIG. 5: The calculated ATR spectrum of a defective multi- 
layer structures for different width of the defect layer dd = 
0.1 di (solid line), dd = 0.8 di (dashed line). Here, the thin 
lines show the case for no = — 1 and the thick lines show the 
case for no = 1. The other parameters are /3 = 1.09, n = 3 
and n = 3. 



uum). We see that when the homogeneous medium is a 
LHM (no = —1), the dispersion cure of the SMs settle 
in the near by the lower edge of the photonic bandgap 
and finally disappears by decreasing the distance of the 
defect layer from the interface of the structure. This is 
in contrast to the case where the LHM is replaced by the 
vacuum. For this case the dispersion cure of the SMs 
moves toward the upper edge of the photonic bandgap. 

As stated in Fig. [21 the localized modes have evanes- 
cent nature and due to this evanescent nature, they will 
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not interact directly with an incoming planewave. So, 
they can be excited by the ATR method. This technique 
has previously been invoked for the investigation of var- 
ious types of surface polaritons, e.g., plasmon-polaritons 
in metals [lol [llj . phonon-polaritons in ionic crystals 
[22L HH, exciton-polaritons in semiconductors I241 . l25ll 
and magnon-polaritons in magnetic materials [26l . |27| . 
We consider the ATR geometry shown in Fig. 1 in ref- 
erence [28j. Here, the uniform medium represents a gap 
layer of the width L that separates dielectric and lay- 
ered structure. For an incident angle larger than the 
angle of the total internal reflection, the electromagnetic 
field incident from an optically dense medium (dielectric) 
with refractive index n > no will penetrate the gap as an 
evanescent wave, which can interact with the evanescent 
localized modes. We have calculated the reflectivity of 
the ATR geometry, using classical electromagnetic the- 
ory. A calculated ATR spectrum for the cases of de- 
fectless and defective multilayer structures are shown in 
Fig. 0] for two different incident angle (or two different 
L3 = nsin(&)). Here, the shaded regions show the first 
spectral bandgap of the multilayer structure. In Figs. 
HJa),(b)) there is only a deep in the ATR spectrum in 
the spectral bandgap corresponding to a SM. However, 
in Fig. [4|c) we reveal two sharp deeps in the ATR spec- 
trum corresponding to a DLM and a SM. In Fig. [U(d), 
we can see only a deep in the spectral bandgap related 
to a DLM. Our investigations reveal that the existence 
region for the SMs and the DLMs not only depends on 
the position of the defect layer (77) but also it depends on 
the thickness of defect layer (dg). To show this, we plot- 
ted the reflectivity of the ATR geometry as a function 
of frequency (k = — ) for different thicknesses of defect 
layer with a fix position in Fig. [5l It is seen that for 
a thin defect layer, the multilayer structure behaves like 
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FIG. 7: The energy flow velocity of (a) SMs and (b) DLMs 
for different thicknesses of the defect layer, dd = 0.1 di (solid 
lines), dd = 0.5 d\ (dotted lines) and dd = 0.8 d\ (dashed 
lines). Here, r\ = 3 and no = — 1. 



as a defectless structure (see solid lines in Fig. [5]) [l8j . 
Increasing the thickness of the defect layer dd, will push 
the SM to the lower (upper) edge of the bandgap in the 
case of no = — 1 ( n o — 1) then it will disappear for the 
sufficiently thick defect layer, and the DLM will appear 
within the bandgap. In the other word, for a big dd the 
structure loses the ability to support the SMs. Therefore, 
the modes localize at the defect layer. 
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FIG. 6: The energy flow velocity of (a) SMs and (b) DLMs for 
different positions of the defect layer, rj = 1 (solid lines), rj = 3 
(dashed lines) and rj — 7 (dotted lines). Here, dd = 0.8 di 
and n = — 1. 



Finally, it is interest of us to know the effect of defect 
layer on the energy flow velocity of the localized modes. 
It is well known that the energy flow velocity vg relates 
total energy flow (S) to the energy density integrated 

over z as ve = jjn- We demonstrated the effect of de- 
fect position and thickness on the energy flow velocity 
of localized modes in the defective multilayer structure 
in Figs. [H and [3 respectively. It is seen that in the 
presence of uniform LHM, the energy flow velocity of 
the SMs can be positive or negative (see Figs. [5Ja) and 
[7{a)), which indicates the considered structure has the 
potential to support forward and backward energy flow 
(this is the noticeable advantage of using LHM). On the 
other hand, Figs, [^b) and \T[b) show that the DLMs 
have only positive energy flow velocity. Accordingly, in 
the used structure the presence of LHM can not affect 
the direction of energy follow of the DLM. In addition, 
Figs. [(Ja) and [7fa) reveal that the presence of defect 
layer can not affect considerably the energy flow velocity 
of the SMs, whereas the energy flow velocity of the DLMs 
depends on the position and thickness of defect layer (see 
Figs, fflb) and Gib)). 
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IV. CONCLUSION 

Briefly, we have studied the dispersion characteristics 
of the localized modes in the presence of defect layer in 
the periodic multilayer structures. It was shown that the 
structure can support both the SMs and the DLMs which 
localize at the interface of the structure and the position 
of the defect, respectively. Moreover, we have revealed 
that the existence regions of the SMs and the DLMs de- 
pend on the position and thickness of the defect layer. 
We show that the dispersion curves of the SMs settle 
near by the lower (upper) edge of the photonic bandgap 
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